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^N Abstract. Asymmetric Dark Matter (ADM) models invoke a particle-antiparticle asymme- 

*y^ try, similar to the one observed in the Baryon sector, to account for the Dark Matter (DM) 

^^ abundance. Both asymmetries are usually generated by the same mechanism and generally 

related, thus predicting DM masses around 5 GeV in order to obtain the correct density. 

The main challenge for successful models is to ensure efficient annihilation of the thermally 

produced symmetric component of such a light DM candidate without violating constraints 

psj from collider or direct searches. A common way to overcome this involves a light mediator, 

'^ into which DM can efficiently annihilate and which subsequently decays into Standard Model 

^ particles. Here we explore the scenario where the light mediator decays instead into lighter 

K> degrees of freedom in the dark sector that act as radiation in the early Universe. While this 

^ assumption makes indirect DM searches challenging, it leads to signals of extra radiation at 

BBN and CMB. Under certain conditions, precise measurements of the number of relativistic 

species, such as those expected from the Planck satellite, can provide information on the 

structure of the dark sector. We also discuss the constraints of the interactions between DM 

and Dark Radiation from their imprint in the matter power spectrum. 



Contents 

1 Introduction 1 

2 Dark Radiation 3 

3 Dark Matter-Dark Radiation interactions 10 

4 Conclusions 16 

1 Introduction 

By now we have overwhelming evidence for the presence of an extra non-baryonic Dark Mat- 
ter (DM) component in the Universe from a variety of different independent sources (see e.g., 
Ref. [1]). These include rotation curves of galaxies, gravitational lensing, structure forma- 
tion and global fits of cosmological data such as the cosmic microwave background (CMB) 
anisotropies measured by the WMAP satellite. However, all this evidence comes exclusively 
from gravitational effects and we remain sadly ignorant of the properties characterizing the 
new particle species constituting the DM, such as their masses or the types and strengths 
of their interactions with the Standard Model (SM) particles. Even in the case where DM 
is composed by a single stable species, it is a reasonable expectation that this will be only 
one ingredient of a richer sector with more particles and interactions, given the complexity 
of the SM sector. The model-building possibilities are therefore limitless, given the drought 
of information on the dark sector (DS) we currently suffer. 

In this context, most of the interest over the last years has been focused in obtaining 
a DM candidate within theories that try to address other shortcomings of the SM. Notable 
examples are the axion [2], introduced to solve the strong CP problem, and the popular 
weakly interacting massive particle (WIMP) [3], with a ~ 100 GeV mass, which can be 
easily accommodated in models addressing the electroweak hierarchy problem. Indeed, the 
masses of the extra degrees of freedom required to stabilize the Higgs mass under radiative 
corrections cannot be much further away from the electroweak scale if the fine tuning problem 
is to be addressed. The fact that the relic abundance of a weakly interacting particle with 
~ 100 GeV mass gives the observed DM energy density is an additional bonus for these 
theories and is usually dubbed the "WIMP miracle" . 

Lately, a paradigm shift in the study of DM models is taking place with the goal of 
exploring the phenomenological consequences of, not only the standard WIMP and axion 
DM models, but also those of other plausible dark sectors. After all, while the electroweak 
hierarchy and the strong CP problems are well-motivated theoretical puzzles on their own, 
the only solid experimental evidence we have so far for physics beyond the SM is DM — along 
with neutrino masses and mixings — and it is fully justified to develop SM extensions with 
only the aim of addressing the existence and nature of DM. This is an important avenue to 
pursue in order to prevent our experimental efforts from becoming too focused on the leading 
DM paradigms, possibly missing relevant phenomenological signals. In this context, the old 
idea of asymmetric dark matter (ADM) [4-9] is becoming increasingly popular [10-64]. ADM 
offers a DM paradigm in which the origin and properties of DM are much more closely related 
to those of baryonic matter. This seems appealing, since both abundances are observed to 



be close to each other Qdm ~ 50^. Indeed, if the origin and mass of the DM candidate 
are similar to the baryons, this coincidence is less striking than within the "WIMP miracle" , 
where the production mechanism and masses in the dark and visible sectors are completely 
different. ADM models postulate that the stability of the DM population stems from a new 
conserved quantum number, X. The relic density is then associated to a particle-antiparticle 
asymmetry, in complete analogy to the baryonic sector and baryon number. A common 
origin for both the baryon and DM asymmetries is usually assumed, which typically implies 
similar abundances and, therefore, a constraint on the DM mass close to 5 GeV so as to 
reproduce the correct observed energy density. This precise prediction turns out to be quite 
general if the mechanism linking the DM and baryon asymmetries conserves a combination 
of B — L and X, say Q. In this case, the DM mass turns out to be 5 — T/Qdm, with Qdm the 
Q-charge of the dark matter particle [65] . The main challenge of successful ADM models is to 
provide sufficient annihilation of the thermally produced symmetric component with such a 
light DM candidate without violating collider or direct search constraints [49, 66]. The most 
common solution involves a lighter mediator with ~ 10 — 100 MeV mass [60] in the DS, into 
which DM can efficiently annihilate and which subsequently decays into SM particles. This 
lighter mediator would be the DS analogue of the pion, which leads to efficient annihilation 
of the symmetric component in the baryon sector. 

In this work we will take the analogy between the dark and visible sectors one step 
further and assume that the DS, in addition to the DM and the mediator, contains very light 
degrees of freedom which, like photons and neutrinos in the SM sector, contribute to the 
radiation content of the Universe. We will dub this content "dark radiation" (DR). Since the 
symmetric component of DM and any mediators can now annihilate or decay into DR and 
not into SM particles, the connection between the DS and the SM weakens, alleviating the 
constraints stemming from collider, DM direct and indirect detection experiments^. Also, in 
models where the symmetry related to the DM number X is explicitly violated by a small 
Majorana mass term, the DM interactions with the DR bath could play an important role in 
DMoanti-DM oscillations, which can normally challenge the survival of the asymmetry [57, 
59, 70]. 

This scenario offers two new probes into the structure of the DS, which are comple- 
mentary to direct, indirect and collider DM searches. On one hand, after the SM and DS 
decouple, the thermal symmetric populations of DS particles end up annihilating via the 
light mediator into DR, which becomes heated with respect to the photon thermal bath. 
Thus, constraints on the allowed amount of extra radiation can lead to constraints on the 
DS degrees of freedom, providing valuable information for model building. The amount of 
energy in dark radiation is traditionally expressed in terms of the extra effective number 
of neutrinos AA^eff that can be probed by its effect on the cosmic microwave background 
(CMB) and the outcome of big bang nucleosynthesis (BBN). Interestingly, these estimates 
of AAeff presently show a trend towards a non-zero DR component [71-84] which shall be 
either confirmed or excluded by the great improvement in sensitivity expected in the near 
future from the CMB data of the Planck satellite. On the other hand, since some interaction 
between the ADM population and the DR must exist via the lighter mediator, this interaction 
can be bounded through constraints on the matter power spectrum. Indeed, for significant 

^However, if the mediator itself constitutes the DR, the long range DM-DM interactions implied are strongly 
constrained through structure formation and can be ruled out in many scenarios [60]. Similar constraints 
have been studied in the context of non-ADM models [67-69] , such as scenarios with a hidden sector photon 
mediating DM-DM interactions [67]. 



DM-DR interactions, DM and DR can form a coupled fluid which, in a way analogous to the 
photon-baryon plasma, is not pressureless and can propagate acoustic oscillations. We will 
discuss the constraints that galaxy surveys can set on this DM-DR interactions. 

The present study on the presence of extra relativistic components in the DS that 
interact at some level with the DM population has been inspired by the asymmetric DM 
paradigm, in which complex dark sectors with extra light degrees of freedom are generally 
required. However, we will try to keep our study as model independent as possible and 
hence our results can also apply to different DM models, not necessarily based in a particle- 
antiparticle asymmetry, in which DM interacts with a DR component. Indeed, the constraints 
on AA'efr stemming from BBN and the CMB analysis that we will discuss can also apply to 
other models of DM. In principle this is also true for the bounds on DM-DR interactions 
that we will derive from the matter power spectrum. However, as we will see, the size of 
the interactions required to lead to any observable effect rules out that DM is a thermal 
relic. Indeed, the annihilation of DM to DR would be too large to reproduce the observed 
DM abundance. The ADM paradigm, on the other hand, decouples the DM abundance, 
determined by a particle-antiparticle asymmetry, from its annihilation cross section and could 
thus lead to the signals we will constrain. Indeed, large annihilation cross sections are 
particularly desirable in ADM models so as to efficiently remove the thermal symmetric 
DM component that can otherwise dominate over the asymmetry and spoil its relation to 
the baryon abundance. 

This paper is organized as follows: In Sec. 2, we derive the constraints that present and 
near future measurements of AN^s imply on the degrees of freedom present in the dark sector 
as a function of the decoupling temperature. Section 3 is dedicated to reviewing the bounds 
that can be derived on the interactions between the DM and DR populations from the galaxy 
power spectrum. Finally, in Sec. 4, we summarize our results and give our conclusions. 

2 Dark Radiation 

If the thermal component of DM ends up annihilating into DR, an extra contribution to 
the energy density in relativistic degrees of freedom will be present, pbr- This contribution 
is generally parametrized through AN^s, the number of extra effective neutrino species by 
normalizing the contribution of the extra radiation to that of a neutrino field 

AiVeff = ^^^ , (2.1) 

■^ 8 30 -^ !^ 

where T^/ is the temperature of neutrinos at the specific moment of interest. 

Many different possibilities can be envisioned for the thermal histories of the SM and 
the DS depending on the details of particular DM realizations, the structure of the DS, and 
its interactions with the SM. A common feature of ADM models is that they usually contain 
significantly more structure than the light stable field constituting the DM. Furthermore, 
given the flavour structure observed in the SM, it seems naive to assume that the DS, which 
amounts to a five times larger fraction of the energy content of the Universe, would actually 
be only composed by a single field. In this context, models with flavoured DM are becoming 
increasingly popular [39, 55, 85-93]. In this work, we will adopt an approach as model 
independent as possible. We parametrize the potential complexity of the DS through the 
number of relativistic degrees of freedom g^ + gi present in the DS at the time of decoupling 
from the SM, characterized by a temperature T^- While gi corresponds to the number of 



degrees of freedom present in the field(s) that ultimately constitute the DR, the Qh degrees 
of freedom correspond to relatively heavy degrees of freedom that are going to turn non- 
relativistic and heat the DR with respect to the SM. This general approach implies that the 
results of this section can apply, not only to ADM models by which we were inspired, but 
also to any DS that contains light degrees of freedom which constitute DR apart from the 
DM candidate. 

Two extreme examples of our parametrization are the following: If only the heaviest 
fields in the DS interact with the SM, T^ will be very high and most of the DS degrees of 
freedom will be relativistic at decoupling, allowing arbitrarily high values of gh depending on 
the complexity of the model. The opposite scenario would be realized when the DR fields 
have interactions with the SM which keep them in thermal equilibrium until very late times. 
In this case, gh = and the DR will not be heated again with respect to the photon bath 
and will have a final lower temperature depending on the moment of decoupling. This last 
example can be realized if DR is made up of sterile neutrinos that are mixed with the SM 
ones and decouple almost at the same time. 

The SM and DS share a common temperature until T^ and from that point onwards they 
evolve independently. In this case, the comoving entropies of the two sectors are conserved 
separately. We can use this fact to track the temperature changes in each sector, which we 
need to evaluate AA'^eff- The comoving entropy in a thermalized sector is defined as 

S=—glT\\ (2.2) 

where T is the common temperature of the sector, a is the scale factor, and gl is the effective 
number of entropy degrees of freedom. The latter can be conveniently expressed as a sum 
over species 

gi{T)= Y. 9^f^ + l E 9^f; (2.3) 

«=bosons j=fermions 

where gi is the number of internal degrees of freedom of species i and 

, _ 45 /SX"^ 4 p yy/y^ - 1 4y^ - I 



are functions of the particle mass m-j through the ratio Zi = nii/T^ that changes smoothly 
from /^ = 1 when the particle is relativistic (zj ^ 1) to /^ =0 when it becomes non- 
relativistic (zj ^1). In practice we can consider them as a filter that allows only relativistic 
species to contribute to g^ . For the SM g^ (T) we have used the approximate fitted expression 
of Appendix A of Ref. [94], which takes into account its non-trivial behavior around the QCD 
phase transition (which cannot be reproduced with the simple formulas above). From now 
on, we will not use any subindex for SM quantities {g^,T) and subindex E8 for DS quantities 

At this stage, it is convenient to develop a general formula for the temperature difference 
of two sectors. After decoupling the comoving entropies are constant during the expansion, 
and so it is their ratio, S1/S2 = gl {Tf / g2,*T2 = constant, where all these four quantities 
can depend on time through the temperatures. The constant can be evaluated at the time 
of decoupling, where Ti = T2 = T^, to be S1/S2 = gl ^{Td) / gl 20^d) ■ The ratio of the two 



temperatures at any later time is therefore 

where the g^^ and Tj are understood to be evaluated at the same time. 

We will now apply this formula to compute AA'gff as a function of the DS degrees of 
freedom. It is convenient to distinguish two different cases. 

High DS decoupling temperature, T^ > MeV 

The simplest and probably more realistic case occurs when the decoupling temperature of 
the hidden sector is higher than the neutrino-electron decoupling, i.e. T^ > MeV. In terms 
of the light and heavy DS degrees of freedom introduced before, the relative temperature 
between DR and photons at the CMB epoch is then given by: 

{9h + gi) 9t:CMB\ Og) 



^7 CMB V 91 9i{Td) 
with S'^cMB being the SM effective number of entropy degrees of freedom at the CMB epoch. 
This number receives contributions from photons and neutrinos, but since the latter ones de- 
couple before the electron/positron annihilation they do not share the corresponding entropy 
injection and have a smaller temperature. The neutrino/photon temperature ratio follows 
from Eq. 2.5 by encompassing photons and electrons/positrons in sector 2, Ty/T^ = (4/11)-*^' '^. 
It follows that ff^cMB ~ ^(■'- + in^) — 3.909. It is worth noting that when electrons 
and positrons annihilate, neutrinos are not fully decoupled and electroweak reactions are 
able to pump a bit of energy into the neutrino sector. This makes the standard value of 
-^cflF — 3.046 [95]. This ~ 2% correction it is too small to be significantly detected even 
by the best prospects of Planck, it teaches us however that incomplete thermalization or 
decoupling of species in the DS can lead to non- integer values of 9h,9i- 

If we assume that all the DS light degrees of freedom have the same final temperature, 
the energy density is /jdr = '^'^9iT^/30, and we find 
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It is easy to prove that in this scenario the value of AA'^gg measured through BEN, A A'^cff | ggj^ , 
is equal or smaller than A A'cff | (jj^g . Before getting into it, we recall that the definition of 
AA^effleBN differs from that of AA^eflflcMB- '^^^ time interval where AA^eglggN affects BBN 
spans from the decoupling of the beta reactions that keep protons and neutrons in chemical 
equilibrium at high temperatures (T ~ MeV) to proper BBN times (end of the deuterium 
bottleneck, T ~ 70 keV). As we already mentioned, between these two boundaries, the 
electron/positron annihilation heats photons but not neutrinos, changing the neutrino to 
photon temperature ratio and therefore the definition of AA'^eff- It is customary to define 
AA^effleBN evaluated at the highest temperature T ~ MeV and this is the definition we use 
in this work. 

At this epoch photons, electrons and neutrinos have the same temperature and 5*33^ ~ 
10.75 and following the same steps as before we obtain 
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where gf^^ is the number of relativistic degrees of freedom in the DS at the time of BBN, 
which may take on values between gi and gi + gh- If gf = gu then AA^eff has the same 
value for CMB and BBN physics, but if larger, the BBN value is smaller. It is interesting to 
note that by fixing gi +gh and T^ we are fixing the comoving entropy of the DS, but the final 
number of degrees of freedom gi (or gf for BBN) still has an impact on AiVefj . Actually, 
the dependence AA^eg oc {gi)~^'^ is quite easy to understand. The reason is that, for a 
fixed entropy density s oc gT^ , the energy density p oc gT^ oc s^'^/(/^'^ is larger for larger 
temperatures and small number of degrees of freedom. This implies that, for two sectors with 
the same number of degrees of freedom at decoupling, that one with the smallest (non-zero) 
number of light species will be more noticeable to AiVcg probes. 

Finally, note that in this scenario we can be sensitive to gf — gi by measuring and 
comparing AiVcfj|ggp^ and AA^cflflcMB- Unfortunately, the foreseeable accuracy of AA^cfrleBN 
makes this comparison challenging. 

Low DS decoupling temperature, T^ < MeV 

In this case, the DS has to couple either to the electromagnetic plasma (electrons, baryons 
and photons) or to neutrinos, which at these low temperatures are decoupled from each other. 
Here we do not consider the case that the DS can mediate interactions between neutrinos 
and the electromagnetic sector. Having a DS in thermal contact with the SM below the 
MeV requires very strong interactions between the two, so the reader should keep in mind 
that many of the models that we can constrain in this section can actually be excluded 
by laboratory experiments or astrophysical arguments. In this case, the DS is still coupled 
during BBN so computing AA'efflBBN amounts to just counting the DS degrees of freedom 
present at T ~ MeV and normalize them to the neutrino energy density 

4 

A^efflBBN = ^(.9H + gh+ gi), (2.9) 

where gn is the number of degrees of freedom that become non-relativistic between the BBN 
and the DS decoupling. We are forced to introduce this new parameter in order to maintain 
the meaning of gt and gi as in the case discussed above. 

Let us now compute AN^sIq^^- In the case in which the DS remains coupled to 
neutrinos until T^ we have 
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where 21/4 are the neutrino degrees of freedom, and we have again assumed T^i < ra^. When 
the DS decouples from neutrinos, it can still get heated with respect to them by the usual 
factor Tjyi/Ty = {{g^ + gi)/ gi)^'^ ■ In this case we find 

^1 (^ , 4 (gfe + ff;)^/3 \ ( ^y,7/A + gH+gh + gi V"' ._.. 

^'='^'^"^=l,'+7^,0^7^jl 3x7/4 + ,. + ,. J • ^^-''^ 

If, on the other hand, the DS couples preferentially to the electromagnetic sector (ES) 
we have T^g = T^ = Td at decoupling. Neutrinos are decoupled from T ~ MeV, where the 
ES degrees of freedom are 2 + 7/2 for photons and e . Thus at T^ they have a temperature 
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where we have assumed T^ <C rUe ~ MeV so that only photons are present at the DS 
decouphng. Also we have used that g'*^(MeV)= gluiTd). Below T^, also the DS decouples 
and later reduces its degrees of freedom from g^ + gi to gi before the CMB epoch. In this 
period the DS temperature increases with respect to the photon one by a factor T^s/T^ = 
{{gh + gi)/giY ■ Since, in this case, the neutrinos do not have their standard temperature 
ratio to photons, it is more convenient to quote the number of effective neutrino species, and 
not only the extra component (AA'^off = -^eflf — 3) 

AT I ,^11 2 + gh + gi \^'\ A ( II gu + gX '^ ,_„. 

Let us note that this scenario suffers from an additional constraint. The entropy contained 
in the H species is dumped into the ES, between the BEN and CMB epochs. This entropy 
ends up in photons at the CMB epoch, but is absent during BBN and therefore the baryon- 
to-photon ratio ry, the parameter that governs the output of BBN will be different between 
these two epochs, 

??BBN ^ 2 + gh+gi 

rjCMB 2 + gH + gh + gi' 
The relic abundance of deuterium is quite sensitive to the value of ?/bbNi and can be therefore 
used to constraint any mismatch between tjcmb and ??bbn, see for instance [96-98]. This is 
because the current BBN prediction using t/bbn = ''?cmb = (6.19 it 0.15) x 10"^'' [99] already 
agrees with the observations within the uncertainties [99]. A nonzero gn will make t/bbn 
smaller than rycMBi increasing the Deuterium yield for a fixed ??cmb- In our scenarios, we 
also have a nonzero AA'eff , but actually this also predicts an increase in D/H so both effects 
go in the same direction. We have checked that indeed these arguments exclude values other 
than gn = 0. 

Bounds, hints and forecasts 

We can now use present data on AA^gff to constrain the structure of the hidden sector, i.e., 
gH,gh and gi. At present, constraints from CMB and BBN physics are similar but not too 
restrictive. There seems to exist a trend in favor of a non-zero value of AN^s > 0, which 
could be hinting at the kind of more complex dark sectors which we study in this work. 

Let us first consider BBN, which has so far provided the strongest constraints on 
AA'efr [100, 101]. An increased value of A^eff leads to a faster expansion of the Universe 
and, as consequences of this, larger Helium mass fraction Yp (since neutrons have higher 
freeze-out abundance and have less time to decay before BBN), larger Deuterium to proton 
ratio D/H (because D burning reactions are less effective) and smaller yields of more massive 
nuclei like Lithium (because they are produced from D and its products, whose reactions are 
slower). At present, comparisons of the observed primordial abundances with the theoretical 
predictions show a trend towards AN^s > 0, with best fit values AA^fj ~ 0.5 — 0.8 [82, 102— 
104]. Interestingly, both the Deuterium and the Helium measurements show this preference. 
The observed primordial abundance of Lithium is much smaller than the predictions from 
standard cosmology. This is the so-called Lithium problem, which can be somewhat alleviated 
by the presence of extra radiation. However, the cited works conclude that these preferences 
are not significant given the errors, and in particular the systematic uncertainties involved in 
the estimation of the primordial abundances. Notwithstanding the above, one can obtain ro- 
bust upper bounds on AA^eff like AA^cff < 1 [102] and AA^cff < 1.26 [104] (both at 95% C.L.). 



These bounds can be relaxed in non-minimal scenarios, as for instance if neutrinos have a 
nonzero chemical potential (up to AA'^eff < 2.56 at 95 % C.L. according to Ref. [104]). Thus, 
it seems that BBN does not contradict the presence of extra radiation, but cannot be used 
to assess it quantitatively nor to exclude it model independently beyond the AA'gfr < 2.56 
level. As we explain in the following, the situation can be much different for the CMB. 

Many claims for an excess AA'gfr ~ 1 during the CMB release are present in the lit- 
erature [71-79] with different levels of significance depending on the dataset and analysis 
performed in the study. While it has been shown that these effects can be amplified by 
volume effects when analyzing data with Bayesian statistics and the effect seems to be sig- 
nificantly prior dependent for some datasets [105], the preference for non-zero AA^gff does 
persist for prior-independent frequentist analyses [80]. Here we will take the latest results 
from the South Pole Telescope collaboration in combination with WMAP data [77] as a 
reference AA^cff = 0.85 ± 0.62. 

It should be noted that the presence of extra radiation during matter-radiation equality 
to which the CMB is sensitive does not necessarily imply its presence during BBN. For this 
reason, we will regard them as independent constraints into two different epochs of the early 
Universe and not combine them. In this work we will mainly focus on the constraints from 
CMB probes. Indeed, the CMB not only offers a window to a lower temperature to which 
DR must contribute even if it was not present at BBN, but also the forthcoming results from 
the Planck mission will soon provide much more stringent constraints superseding present 
BBN sensitivity. Note that, in order to reach its full potential, Planck data should in any 
case be combined with BBN results on the primordial He abundance Yp [106]. 

In Fig. 1 we show the constraints on g^ as a function of T^ for AA'efr = 0.85 it 0.62 
at CMB [77] and a given gi. Below 1 MeV the results correspond to the (less constrained) 
scenario where the DS remains coupled to neutrinos after BBN and gn = 0. For gn 7^ 0, 
similar results are obtained, but interpreting g^ in the vertical axis as g^ + gn- As can 
be seen from this figure, for very late decoupling, having extra heavy degrees of freedom 
is increasingly disfavored. This is due to the fact that the photon bath will not receive 
significant heating from the SM sector at such low temperatures and thus, any extra heating 
in the DS would lead to a too large contribution to AA^eff- On the other hand, if the SM 
decouples from the DS at a higher temperature, then the relativistic degrees of freedom in the 
SM will be heated, requiring heating also in the DS in order for it to contribute significantly 

to AAeff. 

The Planck satellite mission is expected to measure the effective number of neutrino 
species with an excellent accuracy. We compute as a first step the CMB Fisher matrix to 
obtain forecasts for the Planck satellite [107]. Our fiducial model is a ACDM cosmology with 
five parameters: the physical baryon and CDM densities, uj^ and wdm, the scalar spectral 
index, n^, h (being the Hubble constant Hq = 100 h km Mpc~^s~^) and the dimension- 
less amplitude of the primordial curvature perturbations, Ag (see Tab. 1 for their values). 
Furthermore, we add to the ACDM fiducial cosmology a number of DR degrees of freedom 
parametrized as extra sterile neutrino species AA"eflF = 1,2,3. We assume that the ster- 
ile species have thermal spectra and are not coupled among themselves. For these fiducial 
cosmologies, our Fisher forecast analysis provides the following errors: AN^s = 1 it 0.08, 
AAcff = 2 lb 0.08 and AN^s = 3 it 0.1 at la. We then refine this analysis and perform a 
Markov Chain Monte Carlo simulation of the expected Planck results when the total number 
of neutrinos is 3, 4, 5 or 6, which correspond to the cases AN^s = 0,1,2 and 3 respec- 
tively. For the Monte Carlo scan we obtained good agreement with the Fisher matrix results: 



<g 




Figure 1. The la range for the number of heavy degrees of freedom g^ required to heat the light 
sector in order to account for the presently preferred number of extra effective neutrino species during 
CMB AiVoff = 0.85 ± 0.62 and as a function of the decoupling temperature T^. 



ANcs < 0.08, AA^eff = 1 ± 0.10, AiVeff = 2 ± 0.11 and AA^eff = 3 ± 0.14 at la. Therefore, 
near future data from Planck will definitely be able to settle the issue of DR. If evidence 
for iVjy > 3 still persists after these new accurate CMB measurements, it will be extremely 
interesting to further study interacting scenarios in the DR and DM sectors. 
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Table 1. Values of the parameters in the fiducial models explored in this study. 

For the different scenarios that Planck could find, we show in Fig. 2 the vast improvement 
that the smaller errors would imply in the constraints on g^ and, thus, on the complexity 
of the DS. Notice that, in the extreme scenario in which Planck data would prefer 3 extra 
effective neutrino species, some heating from gfi is required even at decoupling temperatures 
as low as 1 MeV and for gi = 3. This extreme scenario is significantly disfavoured by present 
BBN and CMB results, though. In the opposite limit in which Planck finds no evidence 
for extra radiation and only sets an upper limit on AA'^^fj, no contours at the la level are 
allowed for gi = 3 for any decoupling temperature, as long as the SM is only heated by the 
SM particle content with respect to the dark sector after decoupling. 
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Figure 2. The \a range for the number of heavy degrees of freedom g^ required to heat the light 
sector in order to account for a given number of extra effective neutrino species during CMB and as 
a function of the decoupling temperature Td- The bands correspond to different Planck constraint 
forecasts as described in the text. 



3 Dark Matter-Dark Radiation interactions 

The interactions between the DR and DM components will leave an imprint on the galaxy 
power spectrum, see Refs. [108, 109]. In the presence of these interactions, the dark matter 
fluid is no longer pressureless and therefore the situation will be analogue to that of baryons 
and photons before the recombination era, with a series of damped oscillations similar to 
the baryon acoustic oscillations. In analogy to the baryon case, DM-DR interactions will 
modify the matter power spectrum at scales similar to the size of the Universe when they 
become ineffective in changing the velocities of DM particles, after which DM particles begin 
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to fall into potential wells. We denote this typical length scale by 1/fcj ~ l/afH(af) where 
H = dloga/dt is the expansion rate of the universe and the scale factor at freeze-out aj is 
approximately given by solving 

H{af) = r(a/) = nnu (o-dm-dri') = ?^dr (o-oM-DRf^) (3.1) 

Pdm "1dm 

where T is the rate at which DM velocities are changed by 0(1) amounts, p and n denote 
the energy and number density of dark matter or dark radiation, ttidm is the DM mass and 
(-Bdr) the average energy of DR particles. 

It is convenient to parameterize the cross section as 

((TDM-DR-y) = Qo rriBM , (3.2) 

if it is constant or 

Q2 

(dDM-DR-y) = -^ rriBM , (3.3) 

if proportional to T^ [108]. Here Qo and Q2 are constants with units cm^ MeV^^. These two 
cases are representative of possible DR-DM interactions and help to study how the possible 
temperature dependence of the cross section can affect the constraints. The typical scale of 
the DM-DR oscillations in these cases are 



/10-32cm2MeV-i\'/' 
' 0.5 



/ ±u cm ivie V \ , 1 , ^ 

A: ~ 0.5 /iMpc"^ , (3.4) 

\ Qo J 



10-41 cm2 MeV-i ^ ^^^ 



A: ~ 0.6 ( J /iMpc-i , (3.5) 

where we have assumed ANqq = 3. The dependence of these scales on AA'cfj is however quite 
mild, increasing k as the number of effective neutrino species decreases. Figure 3 (upper 
panel) illustrates this effect for a constant interaction cross section, Qq = 10~^^ cm'^ MeV"^, 
for AA'^eff = 1 and AA'efr = 3 interacting with the DM sector. As a comparison, we show 
as well the shape of the matter power spectrum if these species were non interacting. Note 
that the scale at which the damped oscillations appear is well predicted by the approximated 
expression given by Eq. (3.4). The right panel of Fig. 3 illustrates the analogous but for a 
cross section oc T^ and Q2 = 10"^^ cm^ MeV"^. Again, the suppression scale is very well 
approximated by Eq. (3.5). In the following, we shall exploit the dark matter suppression 
effects to set bounds on the interacting rates Qo and Q2 using galaxy clustering data combined 
with other cosmological datasets. 

We have modified the Boltzmann CAMB code [110] incorporating the interacting sce- 
narios and extracted cosmological parameters from current data using a Monte Carlo Markov 
Chain (MCMC) analysis based on the publicly available MCMC package cosmomc [111]. We 
consider here a flat ACDM universe with AA^cff DR species interacting with the dark matter. 
The scenario is described by a set of cosmological parameters 

{u)b,ujr,M, Ho,ns, As,Qo{Q2)} , (3.6) 

where ujh = il^/i^ and wdm = ^dm^^ are today's ratios of the physical baryon and cold dark 
matter densities to the critical density, Hq is the current value of the Hubble parameter, n^ 
is the scalar spectral index, Ag is the amplitude of the primordial spectrum, and Qo; Q2 
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Figure 3. Upper panel: Matter power spectrum for a ACDM model (thick red curve). The blue 
(dotted) dotted-dashed lines depict the matter power spectrum for AN^s = 1 within a (non) interact- 
ing scenario with constant cross section and Qo — 10^'^^ cm^ MeV""'^. The green (long) short dashed 
lines depict the matter power spectrum for AA'off — 3 within a (non) interaction scenario. The lower 
panel shows the analogous but for an interaction cross section c>c 1/a^ and Q2 — 10~^^ cm^ MeV~^. 

encode the DM-DR interactions. Our basic data set is the seven-year WMAP CMB data 
[112] (temperature and polarization) with the routine for computing the likelihood supplied 
by the WMAP team. We analyze the WMAP data together with the luminous red galaxy 
clustering results from SDSS 11 (Sloan Digital Sky Survey) [113], with a prior on the Hubble 
constant from HST (Hubble Space Telescope) [114], including to these data sets Supernova 
la Union Compilation 2 data [115]. 
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Our main results are reported in Tabs. 2 and 3, where we list the constraints on 
the interaction cross section parameters Qo and Q2 in three possible interaction scenar- 
ios, ANes = 1, 2, 3. Notice that, particularly in the more common case of the cross section 
scaling with T^, these constraints are too mild to be significant for the thermal abundance of 
DM. Thus, while the procedure followed is general, the effects studied are mainly relevant for 
ADM scenarios, in which larger than thermal cross sections are required to efficiently annihi- 
late the symmetric component. The bounds are stronger as the number of effective neutrino 
increases, since the typical scale k at which damped oscillations appear in the matter power 
spectrum decreases with Neff. The maximum k used in the analysis of the matter power 
spectrum is ~ 0.15 hMpc~^, since effects that appear at larger scales cannot be easily con- 
strained. Also, note that the constraints found in this study are milder than those found by 
the authors of [108, 109] since in their case the species interacting with the dark matter fluid 
were the active neutrinos and not extra radiation species. There exists a large degeneracy 
between the number of extra radiation species and the dark matter energy density. One of 
the main effects of AiVgff comes from the change of the epoch of the radiation matter equality, 
and consequently, from the shift of the CMB acoustic peaks. The position of acoustic peaks 
is given by the so-called acoustic scale 6a_, which reads 

Oa = ^4M ' (3.7) 

re{Zrec) 

where rg{zrec) and rs{zrec) are the comoving angular diameter distance to the last scatter- 
ing surface and the sound horizon at the recombination epoch z^ec, respectively. Although 
'I'eizrec) almost remains the same for different values of AA'^eff j fsizrec) becomes smaller when 
AA'efr is increased. Thus the positions of acoustic peaks are shifted to higher multipoles i 
(smaller angular scales) by increasing the value of AN^s- The height of the first acoustic 
CMB peak will also increase as AN^s does. Both effects can be compensated by a larger 
cold dark matter component. Therefore, in our analysis, the cold dark matter component is 
larger than in the absence of dark radiation species, and the effect in the suppression of the 
matter power spectrum shown in Fig. 3 will be less noticeable than in the case in which the 
interacting species are active neutrinos and extra DR species are absent. 

Figures 4 and 5 show the la and 2a contours in the (ag, Qo) ^! (^dm^^i Qo) and 
(Age, Qo) planes, with Qo and Q2 in units of 10~^^ cm^ MeV~^ and 10"'*^ cm^ MeV~^, 
respectively, and the age of the Universe in Gyrs. The contours are shown for the three 
possible interaction scenarios considered here, with one, two and three DR species. Notice 
that the scenarios with AA'eff > 1 are increasingly disfavoured. Indeed, taking the minimum 
of the log- likelihood of the different scenarios, we observe a difference of \og{Lmax{ANQff = 
'i)/Lmax{ANcs = 1)) ~ 2. Although this suggests that the scenario with AA'efr = 1 is more 
favored by data, it is not sufficient to rule out a scenario with AA^eff = 3. For this reason 
we included the three scenarios in the figures. Notice that both the physical dark matter 
energy density and the a^ parameter increase as the number of DR species does, since the 
suppression induced in the matter power spectrum by the presence of the extra radiation 
species and also by their coupling to dark matter could in principle be alleviated by a larger 
amount of clustering dark matter as explained above. There exists also a small degeneracy 
between the interaction cross section and the as parameter. This degeneracy can be easily 
understood in terms of Fig. 3 in which we show that the DM-DR interaction decreases the 



^(Tg is defined as the rms matter density fluctuations in spheres of 8 Mpc. 
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amplitude of dark matter fluctuations at small scales. Finally, the age of the universe in 
models with two or three DR species is significantly smaller than in a typical pure ACDM 
universe, since the age of the universe is inversely proportional to the amount of cold dark 
matter in a given cosmological scenario. 

Finally let us comment on the difference between the scenarios with constant and oc T^ 
DM-DR interactions. It is evident from Figs. 4 and 5 that the isocontours very approximately 
transform onto each other if we associate 
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Q, 



Q2 



10-34 c^2 MeV"^ 10-43 cm2 MeV"^ 



(3.8) 



This is not surprising since the DM-DR decoupling is relatively fast for both cases. The two 
cross sections are indeed similar for values of the scale factor aj ~ 2 x 10^, which corresponds 
to the epoch of matter-radiation equality. This is of course, the scale at which the fluctuations 
in the DM density can start to grow fast, unless impeded by the interactions with the DR, 
and sets the characteristic scale for the DM-DR oscillations that we can constrain. 
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Figure 4. The left, right upper panels and the lower panel show the la and 2a contours in the {as, 
Qo), [fluuh'^, Qo) and (Age, Qq) planes, respectively. The interacting parameter Qq is in units of 



10- 



-34 



cm MeV and the age of the universe is in Gyrs. The red, blue and green contours denote 



the three possible interacting scenarios explored here with one, two and three sterile neutrino species 
in the DR sector. 
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Figure 5. The left, right upper panels and the lower panel show the la and 2a contours in the (cg, 
Q2)i (^dm/j^7 Q2) and (Age, Q2) planes, respectively. The interacting parameter Q2 is in units of 
lO^^'^ cm^ MeV^^ and the age of the universe is in Gyrs. The red, blue and green contours denote 
the three possible interacting scenarios explored here with one, two and three sterile neutrino species 
in the DR sector. 



AiV, 



cff 



C.l. 



(95% c.l.) 



1 Qo < 5.6 Qo < 11.8 

2 Qo< 2.6 Qo < 6.2 

3 Qo < 1-6 Qo < 4.6 



Table 2. Upper limits on Qo (in units of 10 ■^ cm MeV ) in different DR scenarios 



AA^, 



cflf 



'0 c.l.) (95% c.l.) 



1 Q2 < 1.8 Q2 < 3.9 

2 Q2 < 1.0 Q2 < 2.7 

3 Q2 < 0.8 Q2 < 1.9 

43 „^ 2 



Table 3. Upper limits on Q2 (in units of 10 ^'^cm^ MeV ^) in different DR scenarios 
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4 Conclusions 

We have studied the phenomenological imphcations of models of dark matter (DM) in which 
the dark sector (DS) also contains additional lighter fields that interact with the DM compo- 
nent and constitute extra dark radiation (DR). We have been mainly inspired in models of 
asymmetric dark matter (ADM) that generally require complex DS extra light components 
into which the symmetric thermal DM abundance can efficiently annihilate. However, we 
tried to keep our analysis as model independent as possible and the constraints we derived 
can also be applied to other models with DM-DR interactions regardless of the asymmetric 
nature of the DM component. Our description involves the splitting of the DS into light and 
heavy degrees of freedom {g^ and g^, respectively). While gi correspond to the DS degrees of 
freedom that will ultimately constitute the DR component, gh parametrizes the DS degrees 
of freedom that were relativistic at the temperature T^, where the DS decoupled from the 
Standard Model, but that became non relativistic and heated the DR at a later time. 

The focus of our study has been put on how cosmological probes on the amount of 
radiation, such as measurements of the cosmic microwave background (CMB), can constrain 
the model. This is particularly interesting given the current preference for additional degrees 
of freedom in radiation displayed by the CMB, as well as the major improvement expected 
in this measurement with the forthcoming release of the Planck results. This provides a 
probe of the DS which is very complementary to more conventional searches, such as direct 
and indirect detection experiments, which would be challenging or absent in this type of 
models. We have found that the DS composition is, at present, a relatively open question 
with up to ~ 20 extra heavy degrees of freedom gh for T^ > 10 GeV. This number is 
significantly suppressed once T^ falls below ~ 0.1 GeV due to the strong reduction of SM 
degrees of freedom that heat the photon bath with respect to DR. Studying the impact of 
the forecasted Planck sensitivity we find that, if the hint for non-zero number of extra light 
degrees of freedom persists, the composition of the DS will be very constrained and this 
could provide important input for DM model building. Furthermore, if the Planck results 
simply put an upper limit, this limit will be strong enough to severely constrain the possible 
existence of light degrees of freedom in this type of models. 

In addition to the effects in the very early Universe, we have also studied the possible 
impact of the interaction between DR and DM on structure formation. If this interaction is 
strong enough, the two would couple allowing the propagation of pressure waves analogous 
to the baryon acoustic oscillations in the baryon/photon plasma and therefore influence the 
galaxy power spectrum. We have studied the bounds on the interactions within the DS (in the 
form of DM-DR scattering) and its correlation with the number of light degrees of freedom. 
While the applicability of this bound is also independent of the asymmetric nature of DM, the 
size of the interactions required for observable effects implies too strong an annihilation of the 
thermal DM component to explain in this way the observed abundance. Thus, this analysis is 
mainly interesting for ADM scenarios in which large annihilation cross sections are desirable 
while the DM abundance is instead controlled by the particle-antiparticle asymmetry, as in 
the baryon sector. We have studied two possible forms of cross sections within the dark 
sector, constant and T'^ dependent. Using Markov Chain Monte Carlo methods, we have 
derived upper bounds on the strength of such interactions, which are summarized in Tabs. 2 
and 3, respectively. We have also seen that cosmological parameters, such as the age of the 
Universe and f^DM^^i have a significant dependence on the number of interacting particles 
that constitute DR in these scatterings. For example, the best fit for the age of the Universe 
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is generally lower in these scenarios as compared to a pure ACDM model. 
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